Tunneling gap of laterally separated quantum Hall systems 
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We use a method of matched asymptotics to determine the energy gap of two counter-propagating, 
strongly interacting, quantum Hall edge states. The microscopic edge state dispersion and Coulomb 
interactions are used to precisely constrain the short-distance behavior of an integrable field theory, 
which then determines the low energy spectrum. We discuss the relationship of our results to the 
tunneling measurements of Kang et al, Nature 403, 59 (2000). 
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For many solid-state materials, modern density- 
functional theory allows accurate computation of elec- 
tronic properties ab initio. This success has not been 
extended to correlated electron systems like transition 
metal oxides, including the cuprate superconductors, and 
low-dimensional semiconductor heterostructures. In such 
systems, the reduced kinetic energy bandwidth enhances 
the importance of the correlations, and in some cases 
this leads to new quantum phases of matter. Quantum 
field-theoretic methods can often lead to a good under- 
standing of the correlated phases, but their quantitative 
predictive power is limited: different energy scales can 
be related to each other, but at least one has to be mea- 
sured experimentally. In this paper we shall describe 
one strongly correlated system in which it is possible to 
make a precise quantitative connection between the mi- 
croscopic Hamiltonian and the electronic spectrum. 

We consider the recent experiments of Kang et 
al. 0| which studied tunneling between the counter- 
propagating edge states of two laterally separated 
quantum Hall states in two-dimensional electron gases 
(2DEG) in GaAs quantum wells (see Fig. |). The 2DEGs 
are separated by a smooth barrier so the elections in 
the edge states move ballistically along the barrier, while 
mixing via a weak matrix element for tunneling under the 
barrier, and strongly interacting via the Coulomb interac- 
tion. These coupled edge states constitute the correlated 
electron system of interest here. As has been argued by 
Kang et al. jj] , and more completely by Mitra and Girvin 
H , the magnetic field and bias voltage dependence of the 
tunneling conductance act as sensitive probes of this sys- 
tem. In particular, the field interval over which a zero- 
bias peak is present determines the energy gap above the 
quantum ground state. 

This paper will make a microscopic computation of 
the energy gap of the coupled edge state system. We 
use the method of matched asymptotics: we combine 
the results of two separate calculations, valid on differ- 
ent length scales, over their intermediate regime of over- 
lapping validity. It is the weak bare tunneling matrix 
element which leads to a clear-cut separation of scales 
crucial to our approach. At short scales, we will use a 




13 um *~ 88A 



14 um 



212 
meV 









_JkL ^ 




_ £ kR 









FIG. 1. The experiment of Kang et al. Jl): two 2DEGs are 
laterally separated by an atomically precise tunneling barrier, 
with a height of 212 meV The differential conductance 
is measured across the junction, with peaks occuring at zero 
bias for filling factors v ~ 1.2-1.5 and v w 2.2-2.5. The edge 
state dispersions are sketched for a barrier of infinite height; 
a gap appears for non-zero tunneling. 



perturbative, microscopic calculation to account for the 
Coulomb interactions between the edge states (in the 
spirit of density-functional theory). At larger scales, we 
will use an integrable quantum field theory to determine 
the energy gap induced by the tunneling. However, the 
results of the quantum field theory depend sensitively on 
the nature of the cutoff regulating its short distance be- 
havior: we will show how this is precisely specified by 
the long distance limit of the first calculation. The same 
quantum field theory has also been used in other recent 
theoretical works ^,|); however, these works have used 
an ad-hoc specification of its cutoff. Our results show 
that there is no ambiguity in regulating the field theory, 
and that there is a precise matching procedure which 
connects it to the microscopic band structure. 

We consider the Hamiltonian H = Hq + H mt + H t . 
For an infinitely high barrier the one-dimensional motion 
of the electrons on the two edge states (a = L,R) is 



1 



represented by 
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where cl is a creation operator for a fermion in an or- 
bital with one-dimensional momentum, k, parallel to the 
barrier. The dispersion was determined by a solution 
of the Schrodinger equation in the lowest Landau level, 
assuming complete spin polarization; we will need the 
full functional form of this dispersion, and not just a lin- 
earized Fermi velocity near the Fermi level. For zero bias 
voltage the dispersion satisfies e/cL = e~kR = hu c v a (k), 
where uj c is the cyclotron frequency. With I the mag- 
netic length, and 2a the width of the barrier, v a {k) de- 
pends on k and a only through V a {k) = w(kl + a/l). 
The single-particle wavefunctions in the Landau gauge 
are i() kcr (x,y) = (f> ka (x) e lkv /\L, where (pka( x ) vanishes 
at the barrier; L is the (effectively infinite) length of the 
barrier. 

The Coulomb interaction among the electrons in the 
edge states is 
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where the matrix elements 
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VkVq = j dx lJ dx 2 <Pl+q <J { x l) ( t ) ka{ x lWc{q,Xi - X 2 ) 

^k'-qa'i^^k'a'i^) (3) 

are computed in terms of the wavefunctions, <j)ka{x), of 
the edge states in the direction transverse to the barrier, 
with 
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Here e is the dielectric constant of GaAs (e = 12.6) 
and qTF is the Thomas-Fermi wavevector accounting for 
screening among the edge-state electrons. Finally, we 
write the tunneling between the edge states as 
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In the absence of interactions this operator causes an 
energy gap of A = 2t in the spectrum of H + H t . 

As t is very small compared to the other energy scales, 
it is useful to first consider the physics for t = 0. The 
system Hq + H m t is a canonical example of a gapless 
Tomonaga-Luttinger (TL) liquid. One of its characteri- 
zations is the power-law correlator 



(O(x)O(O)} = [P , [.I for x 



|4/3 2 



oo, t = 0. (7) 



Here j3 2 is a TL liquid exponent, and A({3 2 ) is a dimen- 
sionless function; both quantities depend upon the com- 
plete microscopic details of Ho + Hi nt . The inverse length 
scale, A, can be chosen at our convenience: only the 
combination A(f3 2 )A 1 ~ 213 is determined by (0), and any 
change in our choice for A will lead to a corresponding 
change in A(0 2 ). None of our results will depend upon 
this choice. For the free particle case (Vc = 0), a simple 
computation shows that (3 2 — 1/2 and ^4(1/2) = l/(27r), 
independent of A. For Vc ^ 0, we computed (3 and A(f3 2 ) 
both by perturbation theory in Vc (valid for ft 2 near 1/2). 
This was the time-consuming part of our calculation, and 
in principle, there is no problem in extending this com- 
putations to higher orders. For now, let us regard /3 and 
A((3 2 ) as known quantities (further details of their com- 
putation appear below) and proceed to a discussion of 
the consequences of H t . 

The effects of H t are more easily discussed in a 
bosonized theory of the TL liquid (t]] . This is expressed in 
terms of the continuum, dimensionless, scalar field </?(x) 
with the Hamiltonian (% = 1) 
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where n is the canonically conjugate momentum to (p, 
and vf is the renormalized Fermi velocity. The standard 
bosonization mapping also shows that the operator 
O(x) is proportional to cos(v / 87r/3<p(a;)), and so we con- 
clude that the long-distance properties of the fermionic 
theory H = H + H int + H t are described by the sine- 
Gordon field theory with Hamiltonian 
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This is a renormalizable and integrable quantum field 
theory, and all its properties are (in principle) deter- 
mined in terms of /3 and a after suitable renormalization 
(normal-ordering) of the cosine interaction. Following 
the convention of Zamolodchikov we normalize the 
cosine interactions by 



cos(V87r/3^(x)) cos(V8tt/3^(0))) 
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for either a — > and all x, or in the short-distance limit, 
x — > and all a The theory (JsJ) contains a soliton 
and an antisoliton with mass M, as well as bound states 
("breathers") with masses m„ = 2M sin(7m£/2), n — 
1, 2, . . . , < l/£ , where £ = f3 2 /{l - (3 2 ), and a and M are 
related by 8 
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Before we can apply ( |TT| ) we need to determine a. It is 
tempting to do this by using the standard representation 
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of the fermionic fields in terms of exponentials of tp and tt 
in the expression for 0{x), and thus map H t to the cosine 
term in H s q. However, this is incorrect: this mapping im- 
plicitly assumes that the fermions are free (Vc = 0), and 
does not properly account for the essential, non-universal 
renormalization of the overall scale of 0(x). Instead, the 
proper answer is obtained with the realization that the 
results (0) and (jll]) in fact refer to the same regime of 
x: the long-distance limit of Hq + iJj nt is described by 
the TL liquid behavior in (Q), and this co-incides with 
the short-distance limit of the continuum theory H s q as 



described by the TL liquid behavior in (10). In this man- 
ner we may deduce the proportionality constant between 
O(x) and cos(v / 87r/3(/?(x)), and so obtain 



a = tA((3 2 )A 
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With the knowledge of (3 and A{j3 2 ) the energy gap, A, 
can be obtained from (]ll])-(|l2|), via 



A = 2min(mi,M) . 



(13) 



In particular A = Ao for the non-interacting case, (3 2 — 
1/2, as expected. 

Before we discuss our final numerical results for A, let 
us present a few details of the promised computation of 
(3 and A([3 2 ). An expansion of the TL correlator (^) for 
P 2 near 1/2 yields 
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On the other hand we calculate the correlation function 
of the tunneling operator in the fermionic theory with 
t = to first order in Hint, 
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and extract the asymptotic behavior of g(x), 

g(x) — ci lnx + C2 + 0(x~ 2 ) for x — > oo . (16) 

Comparison with (|l4j) leads to expressions for 1 and 
A'(^) in terms of c\ and c 2 . It remains to determine C\ 
and C2 from g(x). It is straightforward to calculate the 
required correlation function in terms of a two-particle 
Green function 

(O(x)O(O)) = 

£ e ^(Tc fc , ff ,(r 2 )c fcff (r 3 )4 +gCT (r 1 )4-^(r4)>, (17) 

kk' qaa' 

for vanishing imaginary times, with t\ > t 2 > T3 > T4. 
Up to first order in Vc we obtain 
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where F and its arguments are dimensionless, 
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It is difficult to determine the logarithmic asymptotics 
numerically from this integral over an oscillatory func- 
tion. We therefore consider the Mellin transform G(s) of 
g{x), which is defined by 



G(s) 



dx x s 1 g(x) , s < . 



(20) 



The following particular transform is important for our 
purposes, 



dx x s 
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which shows that powers of logarithms transform into 
poles in the Laurent expansion of G(s). In the present 
case we expect that 



G(s) 



Cl c 2 



+ O{s ) for s . 



(22) 



Performing two partial integration w.r.t q in (|18|) leads to 
the function g(jp), a Fourier-like representation of g{x), 



g(.r)= I dp (I - cospx) d ^ 
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on which the integration in ( |20| ) can be performed, with 
the result (primes mean derivatives w.r.t. p) 
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dpp- s g"(p) + O( S ), 
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where 7 « 0.577 is Euler's constant. For p — > it can 
be shown that g"(p) oc V^qJ} / (vpp); a partial integration 
then yields 



G(s) 
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dp(l-s]np)(pg"(p)y + 0(s°), (27) 



from which we obtain, by comparison with (|22j) , and some 
rearrangement , 
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where an integral representation for 7 was used to sub- 
tract the pole at the origin. The integral in ( p9| ) can be 
performed, with the result 



— = lim 



Cl 



dk 

p^°Jo Jo 



50 d 

dk 'ir 

op 



e -(k+k'+p) 

k + k' +p 



F{-Xk,Xk',X(k + k'+p)) 

A=±l 1 



(30) 



which we evaluate numerically. To obtain the matrix ele- 
ments (j|) we calculated the integral in (||) , which is more 
easily evaluated in terms of a Bessel function, (a, b > 0) 



dp cosp 
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and obtained the real-space wavefunctions of the low- 
est Landau level edge states from the solution of the 
Schrodinger equation. This completes the calculation of 
[3 and A'(l/2). 

Finally, the knowledge of c\ and C2 allows determina- 
tion of the energy gap. Eliminating a from (|ll|)"(^2|) 
and expanding in (3 2 



^, we obtain 
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To this order, it is permissible to use the bare Fermi ve- 
locity w8 = u c £w'(a/£). Matching the expansions in Jl4| ) 
and ( pq ) to express the results in terms of C1.2, the final 
result for the energy gap in the presence of interactions, 
A = 2M, is 
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where the microscopic scale A has cancelled. Note that 
only microscopic parameters appear in this equation; ex- 
cept for the tunneling t itself (which enters via the bare 
gap, A = 2t), they all depend only on properties of the 
two 2DEGs without tunneling. It is also obvious that t 
enters into A in a non-perturbative fashion. 

Our results for A are shown in Table |. We find that 
the correlated gap is larger by a factor of 2-3 than the 
bare gap. Compared with the Hartree-Fock results of Mi- 
tra and Girvin [El our results for A are larger by about 
10 % for the experimental barrier width of 88 A and by 
about 50 % for a smaller barrier of 52 A. We conclude 
that the gap is significantly enhanced by Coulomb corre- 
lations especially for strong tunneling. 

The gap can be related to the filling factor range 
v = v* ± 5v where a zero-bias peak is observed in 
the differential conductance ]l],0]; this signal is expected 
when the Fermi energy falls inside the gap. From our 
results its appearance is predicted in a range of bv « 
A/(huj c ) w 0.01, much smaller than the experimentally 
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TABLE I. Numerical results for the correlated tunneling 
gap A [Eq. (^)] and the TL liquid exponent 2/3 2 , obtained for 
a magnetic field of B = 6 T, corresponding to magnetic length 
i = 105 A and cyclotron energy u c — 10.4 meV (based on the 
effective mass m* = 0.067m for GaAs), and several values 
of the barrier width and Thomas- Fermi screening vector otf . 
The barrier width determines the tunneling gap Ao in the 
absence of interactions and was calculated in Ref. 



observed 8v w 0.15. At present the reason for this dif- 
ference is unclear. It may be due to disorder in the sam- 
ple, which has been suggested to further increase the gap 
0. Another puzzle is that for the non- interacting sys- 
tem the two branches of the dispersion are predicted to 
cross above the second Landau level, so that the zero- 
bias conductance peak is expected near v* m 2, but is 
found experimentally at v* 1 (our approach does not 
provide an estimate of v*). In any case it may be worth- 
while to repeat the experiment with a narrower barrier, 
since in this case the larger bare tunneling amplitudes 
and Coulomb matrix elements should make correlations 
dominate over possible disorder effects. 

In conclusion, we have provided a microscopic and 
quantitative computation of an observable in a correlated 
electron system, the energy gap for lateral tunneling be- 
tween two interacting quantum Hall systems, and com- 
pared our result to experiment. Our method of matched 
asymptotics should find applications in other correlated 
systems. 
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